Quantum brachistochrone problem for spin-1 in a magnetic field 
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We study quantum brachistochrone problem for the spin-1 system in a magnetic 
field of a constant absolute value. Such system gives us a possibility to examine 
in detail the statement of papers [A. Carlini et al, Phys. Rev. Lett. 96, 060503 
(2006), D. C. Brody, D. W. Hook, J. Phys. A 39, L167, (2006)] that the state vectors 
realizing the evolution with the minimal time of passage evolve along the subspace 
spanned by the initial and final state vectors. Using explicit example we show the 
existence of quantum brachistochrone with minimal possible time, but the state 
vector of which, during the evolution leaves the subspace spanned by the initial and 
final state vectors. This is the result of the choice of more constrained Hamiltonian 
then assumed in the general quantum brachistochrone problem, but what is worth 
noting, despite that such evolution is more complicated it is still time optimal. This 
might be important for experiment, where general Hamiltonian with the all allowed 
parameters is difficult to implement, but constrained one depending on magnetic 
field can be realized. However for pre-constrained Hamiltonian not all final states 
are accessible. Present result does not contradict general statement of the quantum 
brachistochrone problem, but gives new insight how time optimal passage can be 
realized. 

Recently Carlini et al. [l| generalized the classical brachistochrone problem to the quan- 
tum case. The quantum brachistochrone problem can be formulated in the following way: 



2 



What is the optimal Hamiltonian, under a given set of constraints such that the evolution 
from a given initial state to given final one \ifif) is achieved in the shortest time. Using 
the variational method the authors solved this problem for some specific examples of con- 
strains which lead to fixed distance between the largest and the smallest energy levels of the 
Hamiltonian. In [2| it was shown that quantum brachistochrone problem can be solved more 
directly using the symmetry properties of the quantum state space. That paper was based 
on the idea considered in [3( , where an elementary derivation was provided for passage time 
from the one quantum state into another orthogonal one. 

Later the variational method was extended to allow finding the time-optimal realization 
of a target unitary operation, when the available Hamiltonians are subjected to certain 
constrains dictated either by experimental or theoretical conditions 4]. In [5| the authors 
considered the brachistochrone problem for quantum evolution of mixed states. Very recently 
Bender et al. studied the brachistochrone problem for a PT-symmetric non-Hermitian two- 
dimensional matrix Hamiltonian jf| and showed that among non-Hermitian PT-symmetric 
Hamiltonians satisfying the same energy constraint the time evolution between two fixed 
states can be made arbitrarily small. Such an interesting phenomenon was observed also 
for dissipative systems described by a non-Hermitian Hamiltonian which has a negative 
imaginary part of eigenvalues Q. Discussion on this subject can be found also in 9]. 

Important statement of work (see also [6|) is that finding the minimal time of 

general evolution reduces to finding optimal time evolution for the Hamiltonian acting on 
the two-dimensional subspace spanned by the initial and final state vectors and \ipf). It 
means that optimal evolution which realizes quantum brachistochrone can be written as a 
linear combination of and \ipf) with time dependent coefficients. One of the aims of our 
paper is to examine this statement in detail within the three-dimensional quantum system. 

We consider the brachistochrone problem in the case when optimal Hamiltonian belongs 
to the pre-constrained class of Hamiltonians, less general, with a less then allowed number of 
free parameters that can be used for the problem. Such case is important from the physical 
point of view when an experimentalist has a possibility to change only a few parameters of 
Hamiltonian but not all. As an example of such a scenario we consider a three level system, 
namely, spin-1 in the external magnetic field described by the Hamiltonian of the following 
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form 



H = hu(n ■ s), 



(1) 



where s are dimensionless spin-1 operators, n is the direction of the magnetic field and fko 
is proportional to the strength of the magnetic field. Eigenvalues of this Hamiltonian are 
— fkv, hu), and 0. The difference between the largest and the smallest eigenvalues is fixed 
AE = HAuo = 2hu, what corresponds to fixed absolute value of magnetic field. 

The Hamiltonian (OQ) contains only two free parameters, namely, two angles 9 and <ft which 
set the direction of the magnetic field 



sin 9 cos ( 



n„ 



sin 9 sin 0, n z = cos 9. 



(2) 



Note that general Hamiltonian in three-dimensional space can be represented by a 3 x 3 
Hermitian matrix which contains nine free parameters (eight if we consider su(3) case). We 
consider here the pre-constrained class of Hamiltonians ([I]) with only two free parameters. 

The brachistochrone problem in this restricted case reads: What is the optimum choice of 
the pre-constrained Hamiltonian, namely, what is the optimal direction of the magnetic field 
n at fixed u, such that the evolution from a given initial state to a given final one \ijjf) 
is achieved in the shortest time? Obviously, with such a restriction of possible evolutions 
not all general final states can be reached. This is the price for taking the narrower family of 
Hamiltonians. However, as we show it below the shortest time achieved in optimal evolution 
is the same as in the general setting despite the fact that evolution in our case is more 
complicated (system leaves the subspace spanned on the intial and final states). Let us see 
it in detail. 

The vector of state for spin-1 contains four parameters. We can write 
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\ 



J 



(3) 



where normalization condition is the following: |a| 2 + |6| 2 -|-|c| 2 = 1. Hence four independent 
parameters, for instance \a\, \b\, a, and a' define the quantum state. Therefore, having 
only two parameters in the Hamiltonian to change we cannot reach arbitrary quantum state 
starting from a given initial one, in other words, the evolution defined by the Hamiltonian (JTJ) 
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cannot relate two arbitrary quantum states. In our restricted quantum brachistochrone prob- 
lem we consider only the states which can be connected by the implemented pre-constrained 
evolution. 

The evolution of the state vector can be realized as follows 



-iHt/h 



- icj(n-s)i 



It is convenient to represent the unitary operator of evolution in the form 



1 — (n ■ s) 2 2 sin 2 i(n ■ s) sin cut. 



(4) 



(5) 



In order to prove this let us note that n • s is the operator of the projection of spin-1 on 
the direction n and it has three eigenvalue —1,0, 1 with corresponding eigenvectors | — 1), 
|0), |1) which can play the role of the basis vectors. An arbitrary vector of state can be 
written as linear combination of these vectors. It is enough to prove formula (jSJ) only for 
basis vectors, which are eigenvectors of n ■ s with eigenvalues —1,0, 1. It is easy to verify 
that for A which takes only three values —1, 0, 1 we have 



1 - A)(l + A) + ^A(A + l)e* + ^A(A - l) e -. 



(6) 



Then using (EJ) for the unitary operator of evolution we just obtain (151). 

Let us take the initial vector of state as the eigenvector of s z with eigenvalue —1 





vj 



(7) 



and reachable final state in form given by ([3]). Then using (jllh and representation ([5]) for the 
operator of evolution, and matrix representation for spin in which s z is diagonal, we finally 
find 

/ --- i2 * sin 2 # sin 2 4 ^ 



\/2, 



e v cos 9 sin 9 sin 



2 uA 
2 



75 c 



' sin 9 sin ut 



\ 



'1 + cos 2 



,2 uit 



sii r ^ + i cos 9 sin ut 



(8) 



/ 



The first component gives the necessary condition that \ip{t)) reaches the final state 

Lot 



sin ^sin z 



a . 



(9) 
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From (jSj) it follows that the second component depends on the first one. Substituting sin 2 ^ 
from ([9]) into the second component of (jHJ) we have 

|&| 2 = 2|a|(l-|a|). (10) 

Then the normalization condition yields the realtion 

| c | 2 = l-|a| 2 -|&| 2 = (l-|a|) 2 . (11) 

Thus we cannot reach arbitrary state, but only such ones which have components satisfying 
conditions ( jTOi) and ( jTTl) . In addition note that the phases for the second and the third 
components are not independent but are related according to (jHJ). If all necessary conditions 
are satisfied, then the time of evolution from the initial state to the allowed final one can be 
found from ([9]) 

4 . ( J\a\\ 
t f = - — arcsm — — - , (12) 
1 Aco \ sin 9 J y ' 

where sin 9 > \/\a~\, KAuj = 2hw is the distance between largest and smallest energy levels. 
It is interesting to note that this expression is very similar to corresponding one for spin- 



1/2 (see, for instance, p, Eq. (5)). The difference is that (1121) contains \/\a\ instead of \a\ 
and is two times larger (in |6j a is denoted as b and Au> is denoted as to). 

We obtain the minimal time which just corresponds to the quantum brachistochrone for 
9 = 7r/2, when magnetic field is perpendicular to z-axis 

Wi = -A- arcsin (^/\a\) . (13) 



Ao; 

As an explicit example let us consider the case \a\ = 1, then |fe| = |c| = and the final 

state 

(l\ 





\ / 



(14) 



is the eigenvector of s z with the eigenvalue 1 and is orthogonal to the initial one (j7|). In this 
case we have the solution for the time only when 9 = tt/2. Thus initial state ([7D evolves 
to final one (IT4"|) only when magnetic field is perpendicular to z-axis. For time of evolution 
we have tf — t min = 2k /Au. This time is two times longer than the shortest possible 
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time obtained in 

HQ 

. Note that the state vector describing the evolution in our case is 
not superposition of initial and final states. Therefore, it is not strange that the time of 
evolution is longer than minimal possible one. 

Let us consider the next example with the initial state 



1 



(15) 



Now the evolution is given by the state vector 
/ 

\m) = 



\e ^ (2 cos 9 sin 9 sin 2 *j + i sin 9 sin uot) 
1 - 2 sin 2 9 sin 2 f 



1 - 2 7 



(16) 



where 



z = —=e 1 ^ ( 2 cos 9 sin # sin 2 — i sin 6 1 sin ut 



i(<^— a) 



(17) 



\z\ =27(1 — 7), 7 = sin ^sin — , tana 



cos(a;t/2) 



2 ' cos 9 sm(ut/2) 

Finally the evolution of the state vector can be represented in the form 



WW) = (1-27) 



+ ^27(1- 7) 



V 





e i(4>-a) 



(18) 



/ 



where 7 and a are functions of time as given in (TTTj) . Let us consider the final state 

f-i\ 

\M = ^= (19) 

1 1 1 

which is orthogonal to the initial one ( fl5l) . In order to reach this state we put 7 (i/) = 1/2. 
This condition gives us the time of evolution 

(20) 



t f = - — arcsm , _ ■ . 
Aw \V2sm9j 



Then choosing additionally that <fi = a(tf) we find that \ip(tf)) = \ipf)- 
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For 9 = 7r/2 we obtain the minimal time of evolution 

tmin = Ao? ( } 

It is interesting to note that this time is equal to the minimal possible time which can be 
obtained according to the statement of [l[ [J where the state vector of evolution for minimal 
possible time belongs to the subspace spanned by the initial and final state vectors or, in 
other words, the vector of evolution for minimal possible time is a superposition of initial 
and final states. In our example, as we see from ( fl8l) . the state vector during evolution does 
not stay in the subspace spanned by the initial and final state vectors. Therefore, we can 
conclude that in order to achieve the minimal possible time it is not necessary that during 
the evolution state vector lies all the time in the subspace spanned by the initial and final 
state vectors. A pre-constrained family of Hamiltonians can yield the more complicated 



evolution, but still with the optimal time, w. 
There is no contradiction with the results of 



rat can have practical value for experiment. 
3], the present result is in agreement with 



the fact that to stay on the subspace spanned on the initial and final states one should use 
the full freedom in the general family of Hamiltonians for the system under consideration, 
but what is important, it shows also that staying on the mentioned subspace, is not crucial 
for implementing the shortest time evolution. 
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